Molecular Dynamics and Structural Studies of Zinc Chloroquine Complexes.

Chloroquine (CQ) is a first-choice drug against malaria and autoimmune diseases. It has been co-administered with zinc against SARS-CoV-2 and soon dismissed because of safety issues. The structural features of Zn-CQ complexes and the effect of CQ on zinc distribution in cells are poorly known. In this study, state-of-the-art computations combined with experiments were leveraged to solve the structural determinants of zinc-CQ interactions in solution and the solid state. NMR, ESI-MS, and X-ray absorption and diffraction methods were combined with ab initio molecular dynamics calculations to address the kinetic lability of this complex. Within the physiological pH range, CQ binds Zn2+ through the quinoline ring nitrogen, forming [Zn(CQH)Clx(H2O)3-x](3+)-x (x = 0, 1, 2, and 3) tetrahedral complexes. The Zn(CQH)Cl3 species is stable at neutral pH and at high chloride concentrations typical of the extracellular medium, but metal coordination is lost at a moderately low pH as in the lysosomal lumen. The pentacoordinate complex [Zn(CQH)(H2O)4]3+ may exist in the absence of chloride. This in vitro/in silico approach can be extended to other metal-targeting drugs and bioinorganic systems.

Water-Triggered, Irreversible Conformational Change of SARS-CoV-2 Main Protease on Passing from the Solid State to Aqueous Solution.

The main protease from SARS-CoV-2 is a homodimer. Yet, a recent 0.1-ms-long molecular dynamics simulation performed by D. E. Shaw's research group shows that it readily undergoes a symmetry-breaking event on passing from the solid state to aqueous solution. As a result, the subunits present distinct conformations of the binding pocket. By analyzing this long simulation, we uncover a previously unrecognized role of water molecules in triggering the transition. Interestingly, each subunit presents a different collection of long-lived water molecules. Enhanced sampling simulations performed here, along with machine learning approaches, further establish that the transition to the asymmetric state is essentially irreversible.

A Blueprint for High Affinity SARS-CoV-2 Mpro Inhibitors from Activity-Based Compound Library Screening Guided by Analysis of Protein Dynamics.

The SARS-CoV-2 coronavirus outbreak continues to spread at a rapid rate worldwide. The main protease (Mpro) is an attractive target for anti-COVID-19 agents. Unexpected difficulties have been encountered in the design of specific inhibitors. Here, by analyzing an ensemble of ∼30 000 SARS-CoV-2 Mpro conformations from crystallographic studies and molecular simulations, we show that small structural variations in the binding site dramatically impact ligand binding properties. Hence, traditional druggability indices fail to adequately discriminate between highly and poorly druggable conformations of the binding site. By performing ∼200 virtual screenings of compound libraries on selected protein structures, we redefine the protein's druggability as the consensus chemical space arising from the multiple conformations of the binding site formed upon ligand binding. This procedure revealed a unique SARS-CoV-2 Mpro blueprint that led to a definition of a specific structure-based pharmacophore. The latter explains the poor transferability of potent SARS-CoV Mpro inhibitors to SARS-CoV-2 Mpro, despite the identical sequences of the active sites. Importantly, application of the pharmacophore predicted novel high affinity inhibitors of SARS-CoV-2 Mpro, that were validated by in vitro assays performed here and by a newly solved X-ray crystal structure. These results provide a strong basis for effective rational drug design campaigns against SARS-CoV-2 Mpro and a new computational approach to screen protein targets with malleable binding sites.